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G\ ■ ABSTRACT 
» 

I Primordial nucleosynthesis constrains the properties of light, stable neutrinos. Apart from 

' the well-known limit on the number of neutrino species, there are also bounds on neutrino 

masses and magnetic moments. I discuss also sterile neutrinos and neutrino propagation in 
a primordial magnetic field, such as could be the origin of the observed galactic magnetic 
> ; fields. 
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1. Introduction 

m 

^ ■ The LEP results have made it clear that there are only three stable light 

left-handed neutrinos which couple to Z. Neutrinos might be exactly massless, 
^1 but even if they have masses, we do not know whether they should be Dirac or 
Majorana masses. For Ue there exist several mass measurements^, all yielding 
"^i/c ~ 10 (although the average (m^ ) is negative). The mass limit on has 
^ ! recently been revised upward^ to 500 keV, although it is likely to be improved 
. soon. The current limit on i^r mass is m,y^ < 31 MeV^. If neutrino masses 
are of the Majorana type, the non-observation of neutrinoless double beta decay 
(Z, A) (Z + 2,A) + 2e- implies that^ (m^) = Y.U^i'm^, < 1-2 eV; barring 
accidental cancellations between the mixing matrix elements Uei, the Majorana 
mass of z^e should thus be less than about 1 eV. 

If neutrinos are Dirac particles, in the Standard Model they will have small 
one-loop induced magnetic moments, given by /Xj, = 3.1 x 10~^^ {m^ / eV) fiB where 
fiB is the Bohr magneton. Studies of z^eC" t'eC" and i'^e~ i'^e~ elastic 
scattering can be used to derive the experimental bounds'^ < 1.1 x 10~^/Ub 
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and jiiy^^ < 7.4 x 10~ /is. The best limit on is obtained from D decays 
in beam dump experiments^, yielding fi^^ < 5.4 x 10~^/Ub. Majorana neutrinos 
cannot have diagonal magnetic moments, but if neutrinos mix, they will have 
non-zero transition magnetic moments. 

There number of astrophysical constraints on neutrino properties, 

mainly based on cooling of stars and the supernova SN1987A. Energy would be 
transferred too fast from the inner core of the supernova if left-handed neutrinos 
flip over to right-handed neutrinos. These would then freely stream out^, thus 
affecting the energetics of the supernova. The most recent numerical study^ 
employs a supernova code that includes also neutrino-nucleon scattering and has 
a higher meson density in the core compared with previous estimates. Therefore 
the importance of processes like tt + N ^ N + is enhanced, and one obtains a 
rather stringent upper limit of 3 keV on the neutrino Dirac mass term. 

There is also the well-known cosmological limit on the sum of all stable 
neutrino masses, whether Dirac or Majorana, which can be obtained by requiring 
that the relic neutrino mass density does not exceed the upper limit of the 
density of the universe. One finds 



where Qoiy is the energy density of neutrinos today, in units of the critical density, 
and h is the Hubble parameter in units of 100 kms~^/Mpc. The age estimates 
of the universe imply^ that floh'^ < 0.25 so that Yli'^'^i ~ 23 eV. If there exists a 
cosmological constant, then the limit is somewhat relaxed; for pvac = O.Spcrit one 
obtains ^ mi^. < 35 eV. 

For cosmology the difference between Dirac and Majorana neutrinos is 
significant because that difference may be crucial for primordial nucleosynthesis 
of light elements^. If neutrinos are Dirac particles, then the right handed 
chirality state vji, as well as the left-handed antineutrino z/^,, must exist. This 
means that at very high temperatures, for each neutrino flavour, there were four 
spin-degrees of freedom in equilibrium (in the Standard Model equilibration of 
right-handed neutrinos takes place e.g. through the one-loop induced magnetic 
moment). Later right-handed neutrino interactions, being very weak, decoupled 
and their number densities were diluted by subsequent annihilations. It is then 
essential that this decoupling occurs at high enough temperature so that the 
relic density of right-handed neutrinos at the onset of primordial nucleosynthesis 
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is small enough. The current nucleosynthesis limit on the maximum number of 
extra degrees of freedom, quantified in units of two-component massless fermions, 
is often quoted as^^ 

SN^ ~ 0.3 (2) 

This obviously sets a limit on the strength of interactions which turn left-handed 
neutrinos into right-handed ones. These interactions will necessarily involve 
spin-flip operators such as the mass and the magnetic moment. 

More generally, nucleosynthesis imposes a constraint on the equilibration rate 
of any sterile neutrino, be it left- or right-handed. This turns out to be a very 
useful way to limit oscillations between sterile and active neutrino species. 



2. Primordial nucleosynthesis 

At temperatures T ^ 0{1) MeV neutrons and protons were kept in chemical 
equilibrium through the weak processes ^ pe, en ^ pUe and n ^ pev^- At 
that time the relative number of neutrons was simply given by the equilibrium 
ratio 

where Am = — rup. The neutron-to-proton ratio froze out at about T ~ 0.7 
MeV, after which the neutron population was still reduced by free neutron decay 
until T ~ 0.1 MeV, at which point photons no longer were energetic enough 
to prevent protons and neutrons to combine to form deuterium. Consequently 
reactions like D + D ^ ^He + n, T + p; D + {p, T, ^He) ^He + 7, '^He + p{n) 
and ^He + ^He '^He + 2p helped to build up about 25% ^He and traces 
of ^He, D and heavier elements such as lithium and beryllium. Therefore, 
assuming the neutron life-time is known, the primordial abundance of '^He 
reflects essentially directly the neutron-to-proton ratio at the freeze-out, which 
depends only on the ratio Vn^p{T) / H{T), where H{T) is the Hubble rate. 

The theoretical Helium abundace thus depends mainly on two things: the 
number of degrees of freedom at the freeze-out which appear in the Hubble 
rate, given hy H = {87r^g,{T) /90y /^TyMpi where g,{T) = Y.9b{T) + |(7f(T) 
is the effective number of degrees of freedom; and the number density of electron 
neutrinos, which can affect the neutron-to-proton rate Vn^p{T). In addition. 
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one needs to know the baryon density of the universe, as weU as the neutron 
hfe-time, which recently has been revised^ to 10.26 ± 0.03 min. There is also 
some inherent uncertainty in the nuclear reaction rates, but for ^ife the main 
uncertainty comes from neutron life-time. A fit to data then yields^*^ 

iV.=3.0-0.SM...(^)-l.(^), (.) 

where Yp is the primordial ^ife-abundance, Tn is the neutron life-time in seconds 
and Tjio < 2.8 is the baryon number in units of 10"''^°. 



Regression analyses of He mass fraction against O and N abundances, with Icr limits. Larger 

and smaller circles represent higher and lower weights, respectively, open circles are objects not 
enriched by Wolf-Rayet stars, and a few typical error bars are shown (from ref. 11). 

The primordial abundances of light elements can be deduced from observations 
in various ways. For instance, one measures the relative ^ife-abundance 
in extragalactic ionized hydrogen regions with different metallicities. The 
result is then interpolated to zero metallicity, which implies^^ a primordial 
abundance Yp — 0.228 ± 0.005 (neglecting possible systematic errors). Adopting 



4 



a conservative upper bound Yp < 0.24 in Eq. (4) gives then rise to the hmit 
Eq. (2). Some quite recent observations of certain very low metaUicity objects 
might however change the fit somewhat^^, and there are also some theoretical 
issues like the nucleon mass corrections on the reaction rates^^, which also alter 
the computed value of Yp slightly. It is nevertheless clear that in any case 
primordial nucleosynthesis precludes the appearance of extra degrees of freedom 
at the level of SN^ <^ 1. 

The fit Eq. (4) can in fact be used to impose a simultaneous bound on the 
neutrino number density and the number of extra degrees of freedom. Let us 
denote the relative i^e -abundance by n^^ and Sriiy^ = riy^ — 1. Then, if g is the 
number of (fermionic) degrees of freedom additional to the Standard Model, one 
can show that nucleosynthesis requires that"*^^ 

g - 4.65n^, < 5N^ , (5) 

with 5Njj deduced from observations as in Eq. (2) or Eq. (4). Electron neutrino 
density may change either due to oscillations or heavy particle decays, which 
take place after the decoupling. 

3. Nucleosynthesis constraints on Dirac neutrinos 

The production rate for the 'wrong helicity' neutrino Vj^, and hence their 
mass, must be small enough so that it decouples already before the QCD 
phase transition which occurs at temperatures somewhere between 100 and 
400 MeV. In that case they will not participate in the entropy transfer from 
quark-gluon plasma to particles in equilibrium, and consequently their number 
and energy densities will be diluted below levels that are acceptable for primordial 
nucleosynthesis. Assuming decoupling just above the QCD phase transition, at 
nucleosynthesis there would then remain a right-handed neutrino energy density 
which would be equivalent to about 0.1 neutrino species. 

The 'wrong helicity' neutrino production rate was first estimated by Fuller 
and Malaney^^, who argued that a Dirac neutrino with a lifetime exceeding the 
nucleosynthesis time scale (t ~ 1 s) should have a mass less than about 300 keV. 
A more detailed study of the scattering processes involved has recently been 
carried out in ref. 16. 
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The starting point is that before the QCD phase transition but below, say, 
T ^ 0.5 GeV, the fermions present in the Universe at significant number densities 
were the leptons and u, d, s and c quarks. There are altogether 47 separate 
2—^2 reactions with (i) no 'wrong-helicity' neutrinos in the initial state and 
(ii) with at least one or z/^ in the final state, which all need to be taken 
into consideration. The constraint (i) \s imposed on because each wrong helicity 
neutrino in the initial (final) state introduces an additional small factor m^/|pp 
(m^/|p'P) to the cross section. Here |p| and |p'| are the absolute values of 
the centre-of-mass momenta of the incoming and outgoing particles, respectively. 
Hence processes with more than one 'wrong-helicity' neutrino can be ignored as 
compared to prosesses with only one 'wrong-helicity' neutrino. In addition, there 
are also quark and lepton decays which can produce i^+'s. 

The thermally averaged scattering rate reads^^ 



K=^,^]:J |^||/(i../T)/(^./T).f-).-(p.„.), (6) 



where is the equilibrium number density of z^+'s, f{Ei/T) are the Fermi-Dirac 
distributions of the incoming particles, and j(pi,P2) = \/ {pi ■ P2Y ~ nalrn^/ E1E2 
is a flux-related factor. In Eq. (6) we have neglected the final state Pauli 
blocking, which is an about 10% effect^^. Thermally averaged decay rate is 
simply given by 



pd _ 1 p(l-234) [ r(-p ,rj..mi /7^ 



where the factor rrii/Ei arises from the Lorentz boost of the decay rate. 

Computing all the relevant processes in this approximation, and requiring 
that r^'^ + < at T = Tqcd, one finds^^ for Tqcd ^ (100)200 MeV that the 
mass limits are m,,^ < (1180)740 keV and m,.^ < (720)480 keV. 

A similar line of argument can be used to set a limit on the magnetic 
moment of a light neutrino. The e'^e~ annihilation cross section for right-handed 
neutrinos in photon-mediated scattering has been estimated to be^^ 

<^-^^(i^V'7l + 8^Vl + 2!=IV (8) 



6mg \ 1 — 4m^ / s 
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Demanding decoupling of the magnetic moment induced interactions prior to 
QCD phase transition yields the limit 



More stringent limits ° will however be obtained from red giants and SN1987A. 



Nucleosynthesis limits on Dirac neutrino masses, as a function of QCD phase transition tempera- 
ture. The allowed region is below the curves. (Prom ref. 16). 

The reasoning described above applies to neutrinos with a mass less than 1 
MeV. A heavy (tau) neutrino with a mass in the MeV region would have a 
more pronounced effect on nucleosynthesis than a light neutrino, because during 
the synthesis of the light elements the energy density of the 'right-helicity' 
states of a heavy neutrino would be comparable to or higher than that of a 
massless neutrino. This is because at that epoch these states have already 
decoupled (T^g^ ~ few MeV). This has been shown^^ to lead to an excluded 
region 0.5 MeV < nii,^ < 30 MeV for the tau neutrino mass, provided r^^ > 10^ 
s. (If 1 s < T^^ < 10^ s, the upper bound is somewhat weakened.) This limit, 
with minor modifications, applies also to Majorana neutrinos. 
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If Vt has a very large magnetic moment, it will be kept in equilibrium by 
photon-mediated annihilations which help to decrease the number density. In 
this way Vt could actually escape the nucleosynthesis constraint above. Following 
the relevant Boltzmann equations one finds^^ in the mass range 5 — 35 MeV the 
bound iii,^ < 10~^, with a slight mass-dependence. The bound is essentially 
determined by the upper limit on ^He. 

4. Sterile neutrinos 

If neutrinos have a non-zero mass, they might also mix with each other 
exactly like the quarks. In the early universe the mixing between flavour states is 
not expected to affect nucleosynthesis^"^, because the number densities of different 
flavour states are (to a high accuracy) equal due to thermal equilibrium. If 
the flavour states however mix with a sterile neutrino z/^, then many interesting 
effects arise. The initial number density of a sterile species may be assumed to 
be diluted by e.g. QCD phase transition, but oscillations may help to flU up 
the density back to its thermal level. Moreover, oscillations may deplete the 
population after the decoupling at T ~ 2.3 MeV but above the n/p freeze-out at 
T ~ 0.7 MeV. Both effects affect primordial nucleosynthesis, as is evident from 
Eq. (5). 

The oscillation pattern between sterile and active neutrinos is in the heat bath 
of the early universe modified by the forward scattering of the active species. 
The average effective energy for v^. at 1 MeV < T <S 100 MeV is given by^^ 

V; = V2GFn^{L - AT^/ml,), (10) 

where ^ ~ 55 and L contains terms proportional to lepton and baryon 
asymmetries; if they are initially small enough, they will be dynamically driven 
to zero^^. Oscillation will not be effective if non-forward collisions destroy the 
coherence, which means that oscillations can start only at temperatures close to 
decoupling. 
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The oscillating system is conveniently described by a 2 x 2 one-body density 
matrix = ^Po{I + P ■ cr), and the equations of motion are^^ 



_ =V X P + 1 - — -^z -(D+ — - " 

dt ^ ^' dt ^ dt 

dn„^ 
~dt 



V X P + (1 - P,)—j^z -{D+ _^)(P,x + Py-y), 
=Fo[Ca{l -ntj- {nt^ - ntj - {nt^ - <J], (11) 



Here a ^ P ^ C {e,/i,T}, Ce = 2.31, C,^ = Cr = 0.51 and Fq = 
2.65 X 10~^(T/McV)^s~^. Vq and V,y^ are the effective energies, generalized from 
Eq. (10) to include also changes in rii,^, D the damping rate and the equation of 
motion for Pq is simply dP^/dt = dn^^/dt. The number densities are normalized 
such that n,^^ = 1 corresponds to a single neutrino degree of freedom in chemical 
equilibrium. A complete description of the evolution of the oscillating system 
Eq. (11) requires numerical study which reveals^^ that the large-angle MSW 
solution for the solar neutrino problem^^ is, in the case of ~ mixing, ruled 
out by nucleosynthesis. Similarly, the nucleosynthesis constraint rules out^^ 
1^1^ — fs oscillations as an explanation for the atmospheric neutrino puzzle^^. (For 
a recent investigation of sterile-active oscillation in the early universe, including 
a numerical nucleosynthesis code, see ref. 30.) 

Sterile neutrinos has also been proposed as a candidate for the cold component 
in a cold and hot mixture of dark matter, which has become popular in the 
view of the COBE detection^^ of the anisotropy in the microwave background. 
It has been argued^^ that sterile neutrino with a Majorana mass in the range 
0.1-1.0 keV, slightly mixed with an ordinary neutrino, would provide warm dark 
matter and structure formation with more power on small scales than in hot 
dark matter scenarios. 

Another interesting possibility, first pointed out by Madsen'^"^ , is to have a 
heavy fermion vh which decays into a light sterile fermion and a light boson. 
These may then be assumed to have decoupled prior to QCD phase transition. 
If uh decay while they are still relativistic, solving for the relevant Boltzmann 
equations one finds that"^^ equilibration in the final state will be preceded by 
decays into very cold {p <^ T) bosons. The bose condensation is effective 
provided the decay rate is fast enough, or equivalently, if the decay temperature 
is large enough: > ^JmHmB where is the boson mass. This would then 
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be a natural way to generate a mixture of hot and cold dark matter. Numerical 
studies sliow'^^ that one obtains about 40% cold bosons, but this could still not 
be enough as the favoured cold/hot ratio appears to be"^^ 70/30. 

5. Spin rotation in magnetic fields 

Right handed neutrinos may also be produced by scattering of i^l off a 
primodial magnetic field, if such exists. This provides an interesting connection 
with neutrino properties and primordial nucleosynthesis on one hand, and with 
the observed galactic magnetic fields on the other hand. The galactic magnetic 
fields, which are of the order of few //G, are believed to have arisen from a 
very weak seed field through the so-called galactic dynamo mechanism"^^ . The 
differential rotation and turbulent motion inside a galaxy amplifies the seed field 
exponentially until a saturation point is reached. The field is observed indirectly 
by measuring the syncrotron radiation of the electrons which traverse the galactic 
field, assuming equipartition of magnetic and particle energies (this assumption 
has however recently been subject to some discussion^^). Observationally not 
much is known about the seed field. Some weak bounds on it may be obtained 
by requiring that the growth time must be longer than galactic rotation period. 
The observation^^ of the magnetic field in a spiral galaxy with z=0.395 seems 
to indicate that the dynamo was saturated already some time ago, implying a 
relatively large seed field. Moreover, the magnetic field of the Milky Way changes 
its direction by about 180° between the Sagittarius and Orion arms, which has 
been argued^^ to be an indication for a large seed field of the order of about 
lO"-*^^ G. In computer simulations^^ seed fields of the order of 10~^^ seem to be 
able to produce the oberved field strength (but not the reversals, which however 
has been observed only in the Milky Way). 

It has been argued^^ that a random seed field of the correct size will indeed 
be produced by fiuctuations in the Higgs field at the electroweak phase transition. 
The point is that at the electroweak phase transition all the physical quantities 
should be uncorrelated over distances greater than the horizon distance. This 
means in particular that the Higgs field cannot be gauge rotated to point into 
same direction in group space in every horizon volume^^. Therefore there 
will appear physical Higgs field gradients from which one may construct an 
electromagnetic field F^- = -i{V^^Vj - ViV^^) with Vi = 2^sm0w/gdi(l)/\(l)\. This 
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is a random field, frozen to the charges in the primeval plasma, and for any 
Gaussian distribution one finds"^^ 

B_.VM=Bo(^)'^, (12) 

where N is the number of correlation lengths and Bq ~ ~ 10^^ G. This 
yields today and at the intergalactic distances of 100 kpc a root mean square 
field of 4 X 10-^9 G with (B) = 0. 

The direct cosmological consequences of the random field Eq. (12) are 
expected to be minor^^, except possibly for Dirac neutrino spin fiip and hence 
nucleosynthesis. For instance, the magnetic energy density at the horizon scale is 
much smaller than radiation energy density. At the onset of nucleosynthesis, and 
at the horizon scale, the magnitude of the field is about Brms — 1500 G. This is 
well below the nucleosynthesis limit B < 10^^ G on primordial magnetic field, 
which is based both on the effects on the expansion rate and on nuclear reactions 
rates. 

To tackle neutrino helicity change in an external slowly varying random 
magnetic field it is best to use a kinetic equation^^ for the Wigner neutrino spin 
distribution function S(p, x, t). It is also useful to consider the combination 

#x(t)e±'"(*) = MBxit) + Ey{t) ± i{By{t) - E,{t)) 



= %j (^[^-(k)e*-^St(k)e-^- ±i(S,(k)e^''-St(k)e-''-)], 

^ By{t)-EM - 

\Ba,{t) + Ey{t). 



(13) 

corresponding to neutrino propagation along z-axis so that 

H±{t) = fi^^B^it) + Blit) + El{t) + El{t) + 2[Ey{t)B,{t) - By{t)E,{t)]. 

The basic starting point for studying neutrino spin fiip in random magnetic field 
is then the neutrino spin kinetic equation for the 2;-component of the neutrino 
spin S{t) — J (i'^pS(p,t), which takes the form^^ 

dS,{t) 



dt 



+ 2#i(t)e^("-^*) J H^{t)e-'^''-^'^S,it)dt 
+ 2^^(t)e-^("-^*) J H^{t)e'^''-^'^S,{t)dt (14) 
= z(c_i^i(t)e-^("-^*) - C+i^i(t)e^("-^*)), 
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and which depends on the fluctuating magnetic field squared. The constants C±i 
are determined from boundary conditions, and H±, a{t) are given by Eq. (13). 
V = 3.45 X 10-2° (T/MeV) MeV is the eflFective neutrino energy as obtained 
from Eq. (10) with L = 0. 

To make use of the complicated kinetic equation (14) one has to average over 
the fiuctuations. This can be achieved by assuming isotropy and using Wick 
rules for the various averages"^^: 

(5t(k)4.(kO) = {27rf6('\k - k'){BlBj), 

Many of the terms in Eq. (14) actually vanish by virtue of the averaging 
procedure. After some manipulations, one finally finds that the kinetic equation 
reduces to the simple expression^^ 

^^^^ 

where the neutrino spin rotation frequency loq is given by 



ujo = ^J{a-VY + AHl. (17) 

Solving Eq. (16) with the appropriate boundary conditions one immediately 
obtains 

l + SJt) 4^2 
P^.^^n = ^ = — ^ sin2(a;ot/2), (18) 

where the frequency now reads 



uo^^V^ + ml) + \L-\ (19) 

Here L is a measure of the magnetic field energy density inhomogeneity, and it is 
determined by the ratio 

A small fluctuation length of the magnetic fleld thus effectively damps the 
neutrino to spin flip probability. Unfortunately, it is not clear what a realistic 
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coherence length of the primordial field should be. Let us however assume that 
above the QCD phase transition the field is coherent over length scales of the 
order of the (electron) neutrino scattering length Lw = (4.0G|.T^)~^ ~ 3x 10~^Ih 
at T = 200 MeV {Ih is the horizon length). We may then average Pl^r over 
Lw to obtain the production rate. One finds that right-handed neutrinos would 
be in equilibrium at nucleosynthesis unless 

r,,^,, = {Pl^r)L^ < H{T = Tqcd). (21) 

In essence, the spin content of the thermal neutrino bath is thus tested by 
ordinary weak collisions. From Eq. (18) one obtains the bound 

//^S(T = 200 MeV, Lw) < 7 x IO^/xbG. (22) 

Assuming the scale dependence of the coherent field is known, we can find the 
magnitude of the galactic seed field at Tqcd- With the scaling Eq. (12) one would 
then obtain the limit fx^ < 2 x 10~^°/iB; with S ~ 1/7V the limit would be^^ 
fJ'v ^ 6.5 X 10~^'^ fj,BG/Bseed{Tnow) — 2 X 10"^^//^. (A similar constraint applies 
also for transition magnetic moments). In the Standard Model this argument 
can also be turned the other way round to yield a limit on the magnetic field 
strength at QCD phase transition, provided galactic dynamo is the explanation 
for the observed galactic magnetic fields: 

B(T = 200 MeV)<^^;^. (23) 

Again depending on the scaling of the magnetic field to intergalactic distances, 
Eq. (23) provides a limit on the sum of all neutrino masses^^. 
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